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The self- consistent solutions of a nonlinear Ginzburg-Landau equations, which de- 
scribe the behavior of a superconducting plate of thickness 2D in a magnetic field H 
parallel to its surface (provided that there are no vortices inside the plate), are studied. 
We distinguish two types of superconductors according to the behavior of their magne- 
tization M{H) in an increasing field. The magnetization can vanish either by a first 
order phase transition (type-I superconductors), or by a second order (type-II). The 
boundary Si-n, which separates two regions (I and II) on the plane of variables (D, x.), 
is found. The boundary (^{p,x) of the region, where the hysteresis in a decreasing 
field is possible (for superconductors of both type), is also calculated. The metastable 
d-states, which are responsible for the hysteresis in type-II superconductors, are de- 
scribed. The region of parameters (D,>c) for type-I superconductors is found, where 
the supercooled normal metal (before passing to a superconducting Meissner state ) goes 
over into a metastable precursor state (p-). In the limit yc l/V^ and D ^ X (where 
A is the London penetration depth) the self- consistent p-solution coincides with the an- 
alytic solution, found from the degenerate Bogomolnyi equations. The critical fields 
Hi, H2, Hp, Hr for type-I and type-II superconducting plates are also found. 

1. Introduction 

The macroscopic Ginzburg-Landau theory of superconductivity (GL) [1] is widely used 
to describe the superconductors behavior in the external magnetic field. In particular, in [1] 
the structure of the intermediate region in the vicinity of superconducting (s-) and normal 
(n-) semi-infinite metallic phases, brought into contact in magnetic field, was studied. It 
was found that the free energy ans of the interface between n- and s-phases vanishes, when 
GL-parameter x equals to xq = 1/ a/2. On this basis all superconductors are usually divided 
in two groops [1]: the superconductors with o"„s > < >fo) belong to first groop, the 
superconductors with ans < belong to second groop. (The value o"„s < indicates that 
s-phase is unstable in respect to the appearance of vortices in a bulk superconductor [2].) 

However, one can use different criterion to divide superconductors in two groops, namely, 
according to the shape of the magnetization curve, using the formula B = H + AnM, where 
B is the mean field value in the specimen, H is the external field, M is the magnetization. 
Such criterion was used in [3] for a cylinder of radius R (assuming no vortices inside the 
cylinder). On the base of the self-consistent solutions of GL-equations the critical parameter 
Xc(-R), dividing two types of the dependence M{H) in the icreasing field, was found. It was 
shown that for h < >Cc (in type-I superconductors) the cylinder magnetization vanishes in a 
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jump (by a first order phase transition to n-state). For x > (in type-II superconductors) 
the magnetization vanishes gradually, by a second order phase transition. (The difference in 
critical parameters xq and Hc{R) is due to the different geometries and different criteria for 
sorting the solutions, assumed in [1,2] and [3].) It was also shown in [3] that the cylinder 
behavior in a magnetic field is by far nontrivial: there exist several states (the solution 
branches), the transitions between these states, the hysteresis phenomena, etc. It was noted 
that the properties of the vortex-free state should be taken into account in interpreting the 
experiments with sufficiently thin (mesoscopic) samples. 

In this paper we consider an infinite plate of thickness 2D in a field H parallel to its 
surface, assuming no vortices inside the plate. (Such vortex-free state was previously studied 
numerically in [4,5], but not in sufficient detail.) We show that even in the simplest case 
of a vortex-free plate the solutions depend on parameters (D, x, H) in a very complicated 
way, analogously to the case of a cylinder geometry [3,6,7]. Thus, there exists a critical 
value Hc{D) ^ xq, which defines two types of the magnetization behavior, M{D, x, H). In 
type-I superconductors (with x < Xc) the magnetization of the plate vanishes (if the field 
increases) in a jump (by a first order phase transition into n-state in some field Hi). In type- 
II superconductors (with x > Xc) the magnetization vanishes gradually (by a second order 
phase transition into n-state in a field H2). If the thickness is sufficiently small, ~ A (A 
is the London penetration depth) , the superconducting state is destroyed by a second order 
phase transition at arbitrary x. If x = Xc(-D), the first and second order phase transitions 
become indistinguishable. We found also, that typc-II superconducting plate, which is in 
the vortex- free Meissner state (with the order parameter ^ ~ 1), may pass (if the field 
increases, Fl-regime) into a special (also vortex- free) "edge-suppressed" e-state [5]. The 
order parameter ip in e-state is strongly suppressed in some layer near the plate surface, so 
the magnetic field B in this layer is not screened and practically equals to the external field 
H. In type-I superconductors such e-state docs not form. (The possibility of transitions into 
e-state means that alongside with the usual vortex mechanism [2] there exists the addidional 
"edge" mechanism of the field penetration into a mesoscopic sample.) 

The hysteresis d-states, which appear in type-II superconductors if the field decreases 
(FD-rcgime), and the so-called "precursor" states (p-) in type-I superconductors, which 
describe the onset of s-state (in the field Hp) from the supercooled n-state, are also studied. 
We show, that in the limit D ^ X and x — > xq = l/\/2 the hysteresis (metastable) p-solution 
coincides with the solution of the degenerate Bogomolnyi equations [8] and can be described 
analytically [9]. 

The paper is divided into several Sections. In Sec. 2 the basic equations and boundary 
conditions are written, and necessary notations introduced; in Sec. 3 the critical lines Si-u, 
( and TT, which exist on a plane of parameters {D, x), are calculated and explained; in Sec. 
4 the superconductor behavior in the vicinity of the critical lines 5'i_n and ( is described; 
in Sec. 5 the examples of the space profiles of different solutions are given; in Sec. 6 the 
metastable solutions, which exist in the region of p-statcs and on its vr-boundary are studied; 
in Sec. 7 the connection of the self-consistent solutions of GL-equations with the solutions 
of the degenerate Bogomolnyi equations is discussed; in Sec. 8 the critical fields Hi, H2, Hp 
and Hr, which characterize the behavior of type-I and type-II superconducting plates, are 
found; in Sec. 9 the results are shortly discussed. 
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2. Equations 



In the case of a plate the system of GL-equations [1] can be written in the following 
dimensionless form: 

(fa 

0, (1) 



(2) 



Here, instead of the dimensioned potential A{x'), field B{x'), current js{x') and the co- 
ordinate x', the dimensionless quantities a(x), h{x) and j{x) are introduced, where 



A-^a B-^b b-— i(x) 



27rA 



3s 



I C(j)o 



(3) 



{x = x'/X is the dimensionless co-ordinate, —Dx < x < Dx, Dx — D/X, X — x^, ^ is the 
coherence length, x is the GL-parameter, (pQ = hc/2e is the fiux quantum). 

We will consider the vortcx-frcc state. In this case the order parameter i/j{x) is an even 
real function, and the potential a{x) (and the current j{x)) are the odd functions of co- 
ordinate: a{x) — —a{—x), i.e. a(0) = 0. Gonsequently, the boundary conditions to Eq. (1) 
can be written in the form (assuming < x < Dx): 



\x=0 





' dx 



(4) 



where hx = H/Hx, Hx = 0o/(27rA2). 

As to the Eq. (2), we will take the usual boundary condition at the external surface 
[1]: difj / dx\x=Dx — 0. The order parameter at the center is maximal, thus, the boundary 
conditions to Eq.(2) are: 



dx 



= 0, 



x=0 



dx 



= 0. 



(5) 



x=D. 



The magnetic moment (or magnetization) of the plate, related to the unity volume, is 
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where Bav is the mean field value inside the superconductor, S is the plate cross-section in 
the (x, y)-plane. In normalization (3), denoting b — Bav/Hx, Mx — M/Hx, one finds: 



AtiMx = b-hx, b 



D: 



-ax, ax = a{Dx), Dx 



D 
J' 



(6) 



The difference of Gibbs free energies of the system in superconducting and normal states, 
AG — Gs — Gn, can be expressed through its magnetic moment: 



HI r 
Stt J 



dv, 



(7) 
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where jF^o corresponds to the superconductor condensation energy, = 0o/(2vrv^A^) is 
thermodinamical critical field. Using (3), one finds from (7) the normalized expression 




9o 



dx 



2^^ + 4 T 



(8) 



The expressions (6)-(8) are used below. 

[Note, that the lengthes ^ and \ — enter the GL-equations on equal footing, each 
one may be chosen as the unit of length; for measuring the field the unit may be either 
H\ = 0o/ (27rA^), or = 4>q/ {2tx^'^) {H^ — >i^Hx). In presenting the numeric results various 
variants of normalization will be used.] 

The solutions were found by using the iteration procedure, described in [10]. In the 
beginning of iterations the trial function was taken as ip{x) ft; 1 in Fl-regime, and %l){x) ~ 0.01 
in FD-regime. The results do not depend on the choice of the concrete numeric algorithm. 

If ^ ^ 1 {B{x) H), the system of equations (1), (2) reduces to a single x-independent 
linear equation 



d'lp/dx^ + (1 - hix')^lj{x) = 



i-D^ <x<D^,D^ = D/i, = H/H^), with ^|;{x) = e"^ w{y), y = 
Weber function w{y) satisfies the confiuent hypergeometric equation w" 
a = 1(1 — l//i^), and can be expressed in a form [11], (2. 44): 

w{y) = CiWi(y) + C2W2(y), 



(9) 

2h^x, where the 
yw' — aw{y) — 0, 



(10) 



^i(z/) = i + E 



a(a + 2) 
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a + 2u-2) 2^ 
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W2{y) = y 



E 



(a + l)(a + 3)---(a + 2i/ 



(21/ + 1)! 

Eq. (9) and the boundary conditions (5) allow to find the critical field of the second order 
phase transition, when ■0(x) — > [12,13]. This critical field does not depend on x. However, 
in the case of first order phase transitions the condition ■?/'(a;) -C 1 is not fulfilled (see below), 
so in a general case to find the critical fields it is necessary to solve the full system of nonlinear 
equations (1), (2), (5). 



3. The state diagram on the plane {D^, x) 



The solutions of Eqs. (l)-(5) depend on the co-ordinate x and three parameters {D\, x, h\). 
In Fig. 1 the state diagram on the plane of parameters {D\, x) is presented. In each point of 
this plane there exist a set of self-consistent solutions for the order parameter iplx; hx) and 
the field h{x; hx). If the representation point {Dx, x) shifts, the corresponding state changes. 
In particular, the magnitude of the order parameter in the origin of co-ordinates changes, 
i/jq = ip{0; hx), and also the mean field value b = hx + AuMx- Such information permits to 
judge how the solution depends on the external field. It is convenient to imagine a peep- 
hole pierced in arbitrary point {Dx, x), what allows to see the dependence of ipo (and of the 
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magnetization, —AnMx = h\ — b) on the field h\ in this point. Studying such dependences, 
one can find on the plane of parameters (Z^a, three critical lines (tt, ^i-h and C), which 
divide this plane into four regions (!„, 1;,, !!„ and lib) with its own characteristic behavior of 
■00 (and —At^Mx) on the field hx. The meaning of these regions is clarified in Fig. 2. 

Fig. 2(a) shows (schematically) the dependence ipo{h\) in region 1^. Evidently, for 
small hx there exists a stable superconducting state with ip^x) ^ 1 {ipQ ^ 1), in which a 
weak external field is almost competely screend and does not penetrate into the bulk of a 
superconductor (the Meissner state, M-). With the field increasing (Fl-regime) the order 
parameter is gradually suppressed, but when the field hx exceeds the critical value hi M- 
solution becomes unstable (small perturbations grow) and passes in a jump {Si) into normal 
state (n-). For hx > hi there is only one stable n-solution with ip = 0. 

Now, searching for the solution in FD-regime with small values of the trial function 
^(x) -C 1, one finds that n-state remains stable at hx < hi down to the restoration field 
hr in Fig. 2(a). In the field interval An — hi — hr two stable states (M- and n-) exist, but 
a supercooled n-state is metastable, because M-state (with ijj ^ 1) has smaller free energy 
(due to the negative condensation energy Qq in (8)). At hx = h,,. n-state looses stability 
(small perturbations grow) with a consequent jump {5r) into M-state. For hx < hr there is 
only one stable M-solution. 

The analogous picture is present in Fig. 2(e), where initial linear grows of the magneti- 
zation {—4:7rMx) corresponds to the Meissner state (6 ~ 0); the transition from M- to n-state 
is acompanicd by a first order jump (Si); there is a hysteresis loop due to the presence of a 
supercooled n-state, and a first order jump (5^) from n- to M-state. 

The characteristic behavior ipo{hx) in region 1^ (Fig. 1) is depicted schematically in 
Fig. 2(b). The picture here is analogous to Fig. 2(a): in Fl-regime there is a jump (6i) 
from M- to n-state; however, a supercooled (metastable) n-state in the field hp passes into 
a special (also metastable) precursor (p-) state, which looses stability in the field h^ with a 
first order jump (5^) to M-state. 

The behavior of the magnetization in this region (Fig. 2(f)) is also characterized by a 
jump (5i); by the presence of a supercooled n-state in the field interval A„ — hi— hp] by the 
presence of p-state (which exists in the field interval Ap = hp — hr), and by a jump (Sr) fromp- 
to M-state in the field hr- The total width of the hysteresis loop is Ap„ = Ap + A„ = hi — hr- 
The width of the interval A^ diminishes in the vicinity to the critical 7r-line (Fig. 1), and 
Ap = on TT-line. Thus, p-states exist only in region I^. 

In region 11^ M-state in Fl-regime (see Fig. 2(c)) becomes unstable in the field hi and 
transforms in a jump (Si) into a new stable "edge-suppressed" state (e-). In this state 
the order parameter is strongly suppressed in some layer near the plate boundary, so the 
magnetic field penetrates this layer practically without screening (see below Fig. 6). If hx 
is further increased, the superconducting e-state is destroyed gradually {ip —>■ 0), with final 
transition to n-state in the field /i2- 

The plate magnetization (Fig. 2(g)) is also characterized by the presence of e-tail and 
by a second order phase transition to n-state in the field /?.2. If now the field is decreased, 
n-state becomes absolutely unstable in the field /i2, where e-statc appears again and for 
fields hx < hi it passes smoothly into a metastable "depressed" (d-) state (see Fig. 6). Such 
d-state (i.e., a supercooled e-state) is a characteristic feature of region 11^, it exists in the 
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field interval A^^ = hi — alongside with M-state and is responsible for a magnetization 
hysteresis (Fig. 2(g)). In region Ila a supercooled n-state is absolutely unstable, in difi^erence 
to region \ (Fig. 2(b,f)) where n-state is metastable in the field interval A„ = /ii — hp. 
If the thickness Dx decreases, remaining in region 11^ (x = const), the jump amplitudes 
5i (and the interval = hi — hj., Fig. 2(c,g)) decrease also, and they vanish (5,,. = Si = 
= 0) on (^-linc (Fig. 1). Below (^-linc (in region II5) the hysteresis is impossible (Fig. 
2(d,h)) and s-state passes into n-state by a second order phase transition. 

The solutions behavior in the vicinity of critical lines tt, Si-u and ( will be studied in more 
details below, here we note the following. As is clear from Figs. 2(e-h), a superconductor 
in the increasing field passes to n-state either by a first order jump, or there is a tail of 
e-states on the magnetization curve and n-state appears by a second order phase transition. 
Acordingly, we will distinguish type-I and type-II superconductors. The boundary between 
two types of superconductors is represented by the curve 5'i_ii (Fig.l). This boundary (or, 
equivalently, the critical value Xc{D\)) depends on the plate thickness, what does not coincide 
with the simple criterion xq = 1/v^ [1], used usually for dividing superconductors into two 
groops. This disagreement (as in the case of a cylinder [3]) is caused by several reasons. 

First, in [1] the case is considered of an infinite superconductor, while we consider a 
superconducting plate of finite thickness. Second, our superconductor boarders a vacuum, 
while in [1] a case is considered of two contacting semi-infinite s- and n-metals. Third, we 
divide two types of superconductors according to the shape of their magnetization curves, 
while in [1] the division is made using different criterion: according to the sign of the surface 
free-energy cr(x) of the interface (with cr = at x = l/\/2 [1]). Thus, the mentioned 
disagreement is due to different settings of the problem. 

4. 5'i_n- and C-boundeiries; G-point 

In this section the results of a self-consistent calculations for magnetization —AT:Mx{h\) 
and ipQ{h\) are given, in a case of a superconducting plate with the parameters Dx, x laying 
near jSi-n-boundary in Fig. 1. 

Fig. 3 represents the case Dx = 7 for three values of x. In Figs. 3(a,b) (x = 0.9, a 
peephole is in region I;,) there arc: M-state, a jump 61 from M- to n-state (in the field hi, 
Fl-regime); a supercooled n-state (A„); a metastable p-state; a jump 6r from p- to M-state 
(in the field hr, FD-regime). 

In Figs. 3(c,d) (x = 0.9193, a peephole is on 5'i_ii-boundary) a supercooled n-state has 
vanished already, but the tail of e-states did not yet appear; this special p-state (or marginal 
^-statc) attains maximum amplitude in FD-regime at the restoration field hr (with a jump 
to M-state). 

In Figs. 3(e,f) (x = 0.95, a peephole is in region 11^) there is a jump {61) from M- to 
e-state, a supercooled n-state is absent, however, a metastable d-state appears which ends 
in a jump (Sr) from d- to M-state in the field hr- 

There is a hysteresis loop on all the curves in Fig. 3; solid lines correspond to Fl-regime, 
dotted lines correspond to FD-regime. 

Fig. 4 illustrates what happens for smaller plate thickness {Dx = 2). The value x = 0.93 
belongs to region I;, in Fig. 1. The value x = 0.953 corresponds to the boundary The 
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value X = 0.98 corresponds to region 11^ (the field interval A^, where a hysteresis rf-state 
exists, diminishes in moving closer to ^-boundary). The value n = 1.03 lies on (^-boundary, 
here the jumps vanish {5i — 5r — A.^ — 0) and the curves become hysteresis-less, having a 
vertical tangent in point io- For x = 1.05 (region llj,) there is no hysteresis, but the inflexion 
point i with finite derivative remains on the curves. If Dx diminishes further, the inflexion 
point lowers and the curves become monotonous, without inflexion. 

It is interesting also to trace what happens in the region of G-point in Fig. 1. The 
critical lines Si^u, C and tt merge in this point {xq ~ 0.915, Dq ~ 1.51) and for Dx < Dq 
there is a single critical curve. Above it (in region !„) the destruction (and nucleation) of 
superconductivity is acompanied by jumps 5i (and 8,). Below it (in region 11^) there is 
a smooth second order phase transition. Thus, for sufficiently small thicknesses all type-I 
superconductors (with x < xq) become, in fact, type-II superconductors. 

To the same conclusion (basing on different considerations) arrived Ginzburg [14], who 
noticed that type-I superconductors (with x <S 1) behave in magnetic field as type-II super- 
conductors. Therefore, G-point may be refered to as the Ginzburg point. [Apart G-point 
there exist the so called tricritical Landau points (L-) [15], where the distinction vanishes be- 
tween the critical fields, which correspond to the supercooled, equilibrium and superheated 
states of a superconductor [16] (i.e., where the hysteresis vanishes). The hysteresys-less 
critical line C, in Fig. 1 consists, in fact, of L-points.j 

Fig. 5(a) shows the dependence oiijjQ and —AtiMx on the field hx in G-point. These curves 
end by a second order phase transition to n-state, having vertical tangent at the transition 
point /i2- Fig. 5b illustates the dependence ipo{hx) in the vicinity of G-point {Dg = 1.51): 
curve 1 corresponds to x = 0.8 (region I^ in Fig. 1); curve 2 corresponds to xg = 0.915 
(G-point); curve 3 corresponds to x = 1.0 (region lib). In region la (curve 1) there exists 
a supercooled n-state (A^J, there are first order jumps {Si and Sr) and the hysteresis loop. 
On curve 2 the derivative dipQ/dhx = oo at the second order phase transition field /i2. The 
hysteresis-less curve 3 with finite derivative at the transition point /i2 corresponds to region 
lib. 

5. Examples of the co-ordinate dependences 

Figs. 1-5 illustrate the solutions behavior on parameters {x, Dx, hx)- In Fig. 6 the 
self-consistent solutions ip{x) and b{x) are depicted as functions of the reduced co-ordinate 
x/Dx, when the representation point crosses the plane of parameters (x, Dx) in Fig. 1 along 
the lines Dx = 7 and Dx — 2. 

Fig. 6a shows the space profile of the order parameter ip^x) for Dx = 7 and several 
valines of x. Curve Me {x = 1.1, region 11^ in Fig. 1) is the Meissner state in the field hi — 
0.9742, which preceeds the jump to e-state {hx — 0.9743). The order parameter of e-curve is 
suppressed near the plate boundary. When the field hx is further increased, the amplitude 
of e-state tends to zero and vanishes finally at /12 = 1.2114. Curve p (x = 0.8, region I5) is a 
precursor state in the field hr = 0.5676, which preceeds the jump to the Meissner state Mp 
{hx = 0.5675). Curve (x = 0.9193, the peephole is on ^i-n-boundary) corresponds to the 
marginal p-state, which attains the maximal amplitude {ipr ~ 0.988) before transforming 
into M^-state (■0o = 0.9999) at hr — 0.6520. The depressed d-state forms in FD-regime from 
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e-state (which exists at h\ = 0.9743, x = 1.1) by gradual transformation of e-solution profile 
and ends (at hr = 0.7797) in a jump to the Meissner M^-state {hx = 0.7796). 
The corresponding profiles of b{x) are depicted in Fig. 6(b). 

The solutions tlj{x) and b{x) for a plate of smaller thickness {Dx — 2) are depicted in Figs. 
6(c,d). Shown are: e-solution at the point of transition from Me- to e-state (x = 0.97, region 
IIq); (i-solution at the point of transition from d- to M^^- {x = 0.97, region Ila); /i-solution 
in the field of transition from /i- to M^- (>«: = 0.953, the boundary ^i-n); p-solution in the 
field of transition from p- to Mp- (x = 0.93, region I;,). The solution is also shown, which 
hes on ^-boundary (x = 1.030, hx — 1.076) with M' = oo at the infiexion point, and the 
solution i {x — 1.05, region 11;,, hx — 1.1) with finite value M' at the infiexion point. 

6. Region of ]>states, 7r-boundary 

In this Section the precursor solutions ip{x) and b{x) are shown at the points (x, Dx) 
belonging to region of Fig. 1. 

The space profiles iIj{x) and b{x) for a plate with Dx = 10 and various x are depicted 
in Figs. 7(a,b). The value x = 0.922 corresponds to ^i-n-boundary (at Dx = 10). At 
this point the jump from M- to n-state happens in the field hi — 0.8511 (Fl-regime, see 
Fig. 3(c)). In FD- regime a supercooled n-state does not exist (it is absolutely unstable), 
but a superconducting p-solution appears, having the amplitude ipQ which grows with the 
field diminishing and reaches the maximum value in the field hr = 0.6529. Such hysteresis 
p-state, which belongs to S'i_n-boundary, will be named as the marginal /x-state. It exists 
in the field interval Ap = hi — hr = 0.1982 simultaneously with M-state and is depicted in 
Fig. 7(a) by curve pL at the field of a jump (hr) to the Meissner state M^. 

If one moves from jS'i_n-boTindary into region I5, a superconducting p-state begins nu- 
cleating from a supercooled n-state at the field hp. The amplitude of p-state reaches the 
maximum at the field hr, after that the jump to M-state occurs. Such hysteresis p-state 
(for X = 0.8, region 1;,) is represented by curve Pr at the point of a jump {hr — 0.5664) to 
Mp-state (not shown); p-state with x = 0.8 exists in the field interval Ap — hp — hr — 0.0744. 

The field interval Ap, where p-state exists, diminishes rapidly with x diminishing. (At 
X = 0.708, region I^, we have hp = 0.5018, hr = 0.5013, Ap = 5 ■ 10^^; at x = 0.707, region 
I(j, we have Ap = 0, with M-state restoring from a supercooled n-state in a first order jump 
without forming p-state preliminary.) The last of p-states existing in region If, corresponds 
to TT-boundary of Fig. 1. Such 7r-state is represented by curve tt in Fig. 7(a) (at the field of 
a jump into corresponding M^r-state, hr = 0.5013). 

The profiles of b{x) for the same states (-Da = 10) are shown in Fig. 7(b). 

If Dx diminishes, the interval of x in region I^ (where the hysteresis p-states exist) di- 
minishes also (Fig. 1) and 7r-boundary curve merges with the curves ^i-n and ( at the same 
point G. There are no p-states in region I^. 

It is interesting to watch how the profiles ip{x) and b{x) of vr-states change, while moving 
along TT-boundary in Fig. 1. This is shown in Fig. 8(a,b) where vr-solutions for Dx = 10 
(x^ = 0.708, hr = 0.5013), Dx = 7 (x^ = 0.708, hr = 0.5014), Dx = 5 (x^ = 0.708, 
hr — 0.5015) and DX — 3 (x — 0.775, hr — 0.6152) are presented. [In distinction to Fig. 
7(b) where b — B/Hx, the fields b^ — B/H^ in Fig. 8(b) are normalized to — x^Hx] 
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The dotted curve in Fig. 8(a) is the solution (W) of the hnear equation (9), normahzed 
to the maximum value of curve 3. It is evident, that the self-consistent solutions ip{x) are 
described at t/jo <til by the Weber functions (10). [Simultaneously, ^ 1 and H fa H^, so 
the linear equation (9) can be used for finding the minimal supercooling field Hr{D), with 
HriD) ^ at D 1]. 

One can see from Fig. 8(a), that when m: l/V^ and Dx ^ 1 the vr-states profiles take 
a characteristic shape of the interface between s- and n- half-spaces [1]. We show below 
that in this special case the metastable (hysteresis) 7r-state coincides with the degenerate 
Bogomolnyi state [8] and can be described analytically [9]. 



7. Connection with the Bogomolnyi equations 



/ 

Jih 



As was shown by Bogomolnyi [8] , at x = 1/ y/2 the GL-equations for the infinite super- 
conductor degenerate and can be reduced to a system of two nonlinear first order differential 
equations, which have the analytic solutions [8,9]. If -0 is a real function only of one Cartesian 
co-ordinate x, the solution is given by the implicit formula [9]: 

'^=J^_ = ±x, b|{x) = l-^^\xl M^) = ^, (11) 

with ■0— >lifx— >— oo, and — > if x — > -|-oo. Point a; = is defined by the condition 
(Pip/dx'^—0 (the infiexion point of ip{x)), i.e., by the equation ■^^ — 1 — lu'^ = with the 
root z^i = 0.451. 

The solid line in Fig. 9 is the solution ip^x) found from the full system of GL-equations 
(TT-solution in Fig. 8(a) at Dx = 10, x — 0.708). The dotted line is the analytic solution 
(11) (— oo < X < +00, X = 1/v^, the infiexion points of both solutions are superimposed). 
Evidently, the self-consistent 7r-solution in the limit Dx ^ 1 coincides with the Bogomolnyi 
solution. The self-consistent field b^{x) (see the solution with Dx = 10 in Fig. 8(b)) also 
matches the formula (11) (6| = 1 — ■0^). [In this connection see [17], where the solutions of 
the Bogomolnyi equations for a single vortex in the infinite superconductor are discussed.] 

In Fig. 7 the profiles of p-solutions ip{x) and b{x) were shown at the point of a jump 
hr from p- to M-state. Fig. 10 illustrates the behavior of p-solutions as function of for 
Dx = 10 and various x. Figs. 10(a,b) demonstrate: (a) - the mean value ip and (b) - the 
free energy Ag (b) in /x-state (x^ = 0.922, the peephole is on the boundary S'i_n). The field 
interval Ap = hp — hr where /x-state exists is Ap — 0.2331 (normalized by H^). In Figs. 
10(c,d) (xp = 0.8, the peephole is inside region I;,) this interval is Ap = 0.1163. In Figs. 
10(e,f) (xtt = 0.708, the peephole is almost on 7r-boundary) Ap = 0.001. 

From Fig. 10 it follows: 1) the degenerate Bogomolnyi solution (B-state, which exists at 
X = l/\/2, Dx oo) is a special case of p-states nucleating in the hysteresis FD-regime from 
a supercooled n-state; 2) B-state exists only in the field hi = hp = hr = 1 (i.e., H — H^) 
when the superconductivity simultaneously originates {ipp ~ 0) and reaches the maximal 
amplitude (V'(O) = 1, but ■0 = 0.5, see Fig. 9); 3) B-state is metastable because M-state 
(-0 = 1) of smaller energy exists as well. 

Notice, that Fig. 10 indicates also to the nonanalyticity of GL-solutions at the degenera- 
tion point (xq = l/\/2, D = oo, h^ = 1). Indeed, at x > xq p-solutions exist, having a form 
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very similar to the degenerate Bogomolnyi solution (with %l) k, 0.5, see Fig. 7). However, 
at X < xo only the absolutely stable M-state remains with a; 1, and absolutely unsta- 
ble n-state with ■0 = 0. In another words, there is a termination point of p-solutions on 
TT-boundary, i.e., nonanalyticity at the point x = xq. 

8. Critical fields (phase diagrams) 

As was mentioned, in an arbitrary point of the state diagram (Fig. 1) the critical fields 
exist {hi and /i2 in Fl-regime, or hp and hr in FD-regime) which are represented schematically 
in Fig. 2. The dependence of the critical fields on the plate thickness may be seen, if one 
makes a mental cut of the plane of states in Fig. 1 along a line x = const. The picture 
seen is presented in Fig. 11 by a number of phase diagrams (in the co-ordinates = H/H^, 

— D /^) ioT plates with different x. 

Dashed line x = 0.3 in Fig. 11(a) corresponds to the critical field hi (Fl-regime); below 
this line lies the region of superconducting M-phase ~ 1), above this line lies the region 
of n-phase {ip = 0). Solid line W corresponds to the critical field hr (FD-regime); above 
this line lies the region of metastable n-phase, below this line lies the region of M-states. In 
region I^ (Fig. 1) a supercooled n-state becomes absolutely unstable in the field hr where 
the jump to M-state occurs. A supercooled n-state exists in the field interval A„ — hi — hr 
where the hysteresis is possible. The interval A„ diminishes with the plate thickness and at 
some the lines hi and hr merge (A„ = 0). This point in Fig. 1 (at x = 0.3) is represented 
by the value Dx = 1.14 {D^ = kD\ = 0.342), which lies on (^-boundary of the hysteresis 
region. For smaller thicknesses (region 11^ in Fig. 1) there exists unique critical field W, in 
which the superconducting state is destroyed (FI) or originates (FD) without hysteresis by 
a second order phase transition. (Notice, that for small x the hysteresis interval A„ seems 
to increase, because the field scale diminishes, = >c^H\.) 

In Fig. 11(a) the critical fields hi and W are also drawn for the values x = 0.5 and 
X = 0.7071 (which lie to the left of the line xq = l/v^ in Fig. 1). Evidently, the critical 
field hi (dashed line) diminishes with x increasing. The fields hr (normalized to H^) are 
represented for all x [12] by a single curve W (which corresponfs to the stability boundary 
of a supercooled n-phase). 

The picture changes in passing to x > xq. Because in region 1;, (Fig. 1) there are 
metastable p-states, we have three critical fields here {hi,hp,hr in Figs. 2(b,f)). These 
fields are represented in Fig. 11(a) (for x = 0.8) by three curves hi,W,hr. The field 
hi (dashed line) corresponds again to the maximum field, at which the jump from M- to 
n-state occurs. At the field hp a supercooled n-state becomes absolutely unstable and a 
superconducing (metastable) p-state of small amplitude originates (this field is presented 
by the same curve W as for x < xq). At the field hr (dotted line) the metastable p-state 
becomes absolutely unstable and a first order transition from p- to M-state occurs. In the 
field interval Ap = hp — hr the metastable p-states exist and the hysteresis is possible. 

The field interval Ap diminishes with the plate thickness, so crossing vr-boimdary (Fig. 
1) we get from region I^ into region I^ where p-states are absent but the supercooled n-state 
exists. (For x = 0.8 7r-boundary is crossed at Dx — 2.7 or — 2.16. The dotted fine hr 
in Fig. 11(a) merges with the curve W at point a.) In region I^ there are only two critical 
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fields {hi and W). Further decreasing D we get from region into region II5 crossing tlie 
liysteresis ^-boundary (for x = 0.8 tliis liappens at = 1.2, point /? in Fig. 11(a), when 
the curves hi and W merge). In region 11^ {D^ < 1.2) no hysteresis is possible and in Fig. 
11(a) remains only one critical curve W which describes the states with t/j ^ 0. 

In Fig. 11(b) the critical fields are depicted for x = 1; 1.2; 2. If x = 1 (the mental cut 
in the plane of Fig. 1 lies inside region 11^) there are three critical fields: hr, hi, h2 (see Fig. 
2(c,g)). The tail of e-states with 1 either vanish (FI) or appears (FD) in the field h2- 
This field is marked in Fig. 11(b) by a letter W (solid curve). In the field hi (dashed fine) 
the jump from M- to e-state occurs (in Fl-regime) and d-state appears (in FD-regime). The 
metastable d-states (and the corresponding hysteresis) exist down to the field hr where the 
jump from d- to M-state occurs. The width of the hysteresis region A^^ = hi — hr depends on 
the plate thickness. The points A,^ = in Fig. 11(b) (or the Landau points, L) correspond 
to the intersection of the line x = const with ^-boundary in Fig.l. (For x = 1 A^ = at 

— 2.0.) In region Ilf, the hysteresis is absent, so here exists only one critical field W 
which describes the reversible second order phase transition. 

Phase diagrams at x = 1.2 and x = 2 are analogous to the case x = 1. The second order 
phase transition curve W (as well as the curve W in Fig. 11(a)) is the same for all x, it 
can be found from the linearized equation (9) and expressed through the Weber functions 
w (10). However, to find the first order phase transition fields {hi and hr) it is necessary to 
solve full system of GL-equations. 

Note, that the interval between L-point (where A^^ = 0) and W^-curve diminishes with x, 
so the curves of Fig. 11(b) transform continously into the curves of Fig. 11(a). 

9. Conclusion 

Note in conclusion, that the vortex-free states, studied above, may be realized in meso- 
scopic samples with characteristic dimension D of several A. With D increasing the uniform 
(one-dimensional) edge-suppressed e-state (as well as d- and p-states) may become unstable 
relative breaking the boundary region into separate vortices, with forming subsequently the 
regular vortex lattice [2]. However, the detailed study of such inhomogeneous states demands 
the solution of partial differential equations what is outside the scope of the present inves- 
tigation. (In this connection see, for instance, Refs. [18,19] where some of such problems 
are considered to explain the experiments [20-24] with thin superconducting discs of various 
form in a perpendicular magnetic field.) Besides, many of the theoretical results happen 
to depend rather weakly on the sample geometry, so the predictions obtained for the plate 
(or the cylinder [3]) on the base of one-dimensional equations, have, probably, more general 
value and may be used in discussing the details of the concrete experiments. 

I am gratefuU to V. L. Ginzburg for the interest in this work and valuable comments, 
and also to A.Yu. Tsvetkov and V. G. Zharkov for discussions. This work was supported 
through the grant RFFI 02-02-16285. 



11 



References 



[I] V.L.Ginzburg, L.D.Landau, Zh.Exp.Teor.Fiz. 10, 1064 (1950). 

[2] A.A.Abrikosov, Fundamentals of the Theory of Metals (North- Holland, Amsterdam, 
1988). 

[3] G.F.Zharkov, Zh.Exp.Teor.Fiz. 122, A^9 (2002). 
[4] P.M.Markus, Rev. of Mod. Phys., 36, 294 (1964). 

[5] A.Yu.Tsvetkov, G.F.Zharkov, V.G.Zharkov, Krat.Soob.Fiz. FIAN, A^2, 42 (2002). 
[6] G.F.Zharkov, V.G.Zharkov, A.Yu.Zvetkov, |cond-mat/00082171 (2000); 

Phys.Rev.B, 61 12293 (2000). 
[7] G.F.Zharkov, JLTP, 128, A^3/4 (2002). 

[8] E.B.Bogomolnyi, Yad.Fiz. 24, 861 (1976) [Sov. J.Nucl.Phys. 24, 449 (1976)]; 

E.B.Bogomolnyi and A.I.Vainstein, zbzd. 23, 1111 (1976) [23, 588 (1976)]. 
[9] I.Luk'yanchuk, Phys.Rev.B, 63, 174504 (2001). 
[10] G.F.Zharkov, V.G.Zharkov, Physica Scripta 57, 664 (1998). 

[II] E.Kamke, " Differentialgleichungen" , part I, Leipzig (1959). 
[12] D.Saint- James, P.deGennes, Phys.Lett. 7, 306 (1963). 

[13] D.Saint- James, Phys.Lett. 15, 13 (1965). 

[14] V.L.Ginzburg, Zh.Exp.Teor.Fiz, 34, 113 (1958) [Sov.Phys. JETP, 34, 78 (1958)]. 

[15] L.D.Landau, Phys.Zs.Sowjetunion 8, 113 (1935); 11, 26 (1937). 

[16] H.J. Fink, D.S.McLachlan and B.Rothberg-Bibby, in Prog, in Low Temp. Phys., 

V. Vllb, p. 435, ed. D. Brewer, North Holland Pub., Amsterdam- N.Y.- Oxford (1978). 
[17] V.Hakim, A.Lemaitre, K.Malhck, Phys.Rev.B, 64, 134512 (2001). 
[18] V.A.Schweigert, F.M.Peeters et.al., Phys.Rev.Lett. 79, 4653 (1997); 

ibid, 83,2409 (1999); Supralatt. and Microstruct. 25, 1195 (1999); 

Phys.Rev.B 59, 6039 (1999); ibid, 62, 9663 (2000); Physica C, 332, 266,426,255 (2000). 
[19] J.J.Palacios, Phys.Rev.B 57, 10 873 (1998); Physica B, 256-258, 610 (1998); 

Phys.Rev.Lett. 83, 2409 (1999); 84,1796 (2000). 
[20] V.V.Moshchalkov et al.. Nature (London) 373, 319 (1995); 408, 833 (2000). 
[21] A.K.Geim et al.. Nature (London) 390, 259 (1997); 396, 144 (1998); 407, 55 (2000); 

Phys.Rev.Lett., 85, 1528 (2000); 86, 1663 (2001). 
[22] D.S.McLachlan, Solid State Commun., 8, 1589, 1595 (1970). 
[23] O.Buisson et al, Phys.Lett.A 150, 36 (1990). 
[24] F.B.Muller-Allinger, A.C.Motta, Phys.Rev.B 59, 8887 (1999). 



12 



Figures captions 



Fig. 1. The state diagram on the plane {Dx, x). Curve 5'i_n is the boundary between 
first and second order phase transitions from s- to n-state in the increasing field (Fl-regime); 
curve 71 is the boundary of metastable p-states in the decreasing field (FD-regime); curve ( is 
the hysteresis boundary, below ^-line the hysteresis is absent. The asymptotics of ^-boundary 
are: = 1.13 for x ^ 0; Da = 2.43 for x > 6. 

Fig. 2. The order parameter t/jQ and magnetization {—AnMx) versus field hx in different 
regions of Fig. 1 (schematically, see the text). 

Fig. 3. The dependences of ipQ and —AtiMx on hx in the vicinity of S'i_n-boundary in 
Fig. 1 [Dx = 7, the values of x are given in the figure.) M is the Meissner state; p is the 
precursor state (region I^); is the marginal p-state (laying on the critical .Si-n-boundary); 
d is the metastable depressed d-state (region Ila). Dotted fines correspond to FD-regime. 

. 4. Tfie dependences of ipQ and —AttMx on hx at Dx = 2 and various x: x = 0.93 
(region I;, in Fig. 1), >«: = 0.953 (/x-state on S'i_n-boundary), >«: = 0.98 (region 11^), x = 1.03 
(^-boundary), x = 1.05 (region 11^); p, /i, d are tfie fiysteresis (metastable) states. On tfie 
curves x — 1.05 tfie points of inflexion i are marked. 

. 5. (a) - Tfie dependences ■0o(^a) and —AnMxihx) in G-point in Fig. 1. (b) - Tfie 
dependences ipo{hx) in tfie vicinity of G-point (Dq — 1.51): 1 - h — 0.8 (region Iq); 2 - 
x= 0.915 (G-point); - x= 1.0 (region 11;,). 

. 6. (a) - Tfie order parameter '^{x) and (b) - tfie field hx{x) [Dx = 7) in different states: 
e - X = 1.1, hi = 0.9743; p - x = 0.8, hr = 0.5676; - x = 0.9193, hr = 0.6520; d - x = 1.1, 
hr — 0.7796. Tfie corresponding M-states (see tfie text) are also sfiown. 

(c) and (d) - Tfie analogous curves for Dx = 2: e - x = 0.97, hi = 1.0172; p - x = 0.93, 
hr = 0.9612; - x = 0.953, = 0.9903; - x = 0.97, hr = 1.072; io - x = 1.030, 
hx = 1.076; i - >t: = 1.05, hx = l-l (see tfie text). 

Fig. 7. Tfie precursor states: (a) - ip{x) and (b) - bx{x) for Dx ~ 10 and different x: 
- (,Si_ii-boundary) x = 0.922, hr = 0.6529; p^ - x = 0.8, hr = 0.5664; tt - x = 0.708, 
hr = 0.5013 (see tfie text). 

Fig. 8. Tfie space profiles of p-states, laying on vr-boundary at tfie field of a jump (hr) 
from p- to M-state (for Dx = 5; 7; 10). Tfie values of x and hr (normalized to Hx) are given 
in tfie text. Tfie dotted fine in (a) is tfie normafized solution W of tfie finear equation (9). 

Fig. 9. Sofid lines are tfie self-consistent solutions for ■0(x) and b{x) (normalized to H^) 
for Dx — 10, K — 0.708. Dasfied line is tfie degenerate Bogomolnyi solution (11). Tfie 
solutions are superimposed at tfie inflexion points tjji = 0.451. 

Fig. 10. Tfie mean value ijj (a) and tfie free energy Ag (b) as functions of h^ for tfie 
states existing in tfie plate of tfiickness Dx = 10 and various x (sfiown in figure). Solid 
lines are - and n-states, dasfied lines are tfie precursor states (//, p, tt) . Points B on vr-curves 
correspond to tfie degenerate Bogomolnyi solution witfi ijj ~ 0.5. 

Fig. 11. Tfie critical fields for different x (sfiown in figure): (a) - type-1 superconductors, 
X < Xc, (b) - type-II superconductors, x > Xc. Dasfied lines are tfie fields hi in wfiicfi M- 
state becomes absolutely unstable (Fl); sofid fine W is tfie field in wfiicfi n-state becomes 
absolutely unstable (FD); dotted lines are tfie fields hr in wfiicfi M-state restores (FD). 
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